Discrete Breathers at the Interface 
between a Diatomic and a Monoatomic Granular Chain 
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In the present work, we develop a systematic examination of the existence, stability and dynamical 
properties of a discrete breather at the interface between a diatomic and a monoatomic granular 
chain. We remarkably find that such an "interface breather" is more robust than its bulk diatomic 
counterpart throughout the gap of the linear spectrum. The latter linear spectral gap needs to 
exist for the breather state to arise and the relevant spectral conditions are discussed. We illustrate 
the minimal excitation conditions under which such an interface breather can be "nucleated" and 
analyze its apparently weak interaction with regular highly nonlinear solitary waveforms. 



INTRODUCTION 



Over the past decade, granular chains have at- 
tracted significant attention [iHl]- This is because 
of the unusual nature that allows the tailoring of 
the response of the chain to be at the essentially lin- 
ear, weakly nonlinear or strongly nonlinear regime 
0. In particular, these chains of solid particles in 
the absence of static precompression are subject to 
purely nonlinear elastic contact interactions between 
their constituents leading to fundamentally novel 
features such as a vanishing sound speed (the so- 
called "sonic vacuum") [l[. This renders the linear 
and weakly nonlinear continuum approaches based 
on the Korteveg-de Vries (KdV) equation Q invalid 
and places granular chains in a strongly nonlinear 
regime of wave dynamics. Such a granular arrange- 
ment supports a new type of nearly compact highly 
tunable solitary waves that have been experimen- 
tally and numerically observed in one-dimensional 
(ID) Hertzian granular systems (sl-flTj. While spher- 
ical bead chains have been mostly studied, similar 
responses can be obtained from cylindrical and ellip- 
tical particles, and fibrous layers fl^, among others. 

More recently, it has been recognized that hetero- 
geneous type settings bear significant additional po- 
tential for interesting nonlinear excitations. In par- 
ticular, in diatomic, triatomic etc. chains, the pos- 
sibility of appropriately modified, nearly compact 
traveling waves was revealed in [l^, and [1^; on 
the other hand, more recently, the ability of such di- 
atomic chains of spheres with Hertzian interactions 
under precompression to sustain discrete breather 
(DB) type excitations was proposed and experimen- 



tally corroborated. The theoretical [21[ and experi- 
mental realization [2^ of such nonlinear vibrational 
structures followed up on earlier works that had ar- 
gued the emergence of nonlinear, defect-induced vi- 
brational modes [2^, even in the absence of pre- 
compression [2J]. In such acoustic vacuum settings, 
breathers have been argued to be metastable [2^ ; it 
should also be mentioned that the linear dynamics 
of such diatomic crystals has been explored in the 
case of welded sphere contacts in [2^ . 

It is at the latter junction of granular chains and 
nonlinear localized vibrations that the present work 
is centered. In particular, we were motivated by 
the recent experimental investigation of 27[ exam- 
ining the propagation of excitations at the interface 
between a monoatomic nonlinear chain and a lin- 
ear medium to consider the possibility of formation 
of vibrational excitations at an interface. In Refs. 
(2lj, and [2^. we identified such modes in the "bulk" 
of the diatomic chain. Specifically, we found two 
classes of solutions: (1) an unstable discrete breather 
that is centered on a heavy particle and character- 
ized by a symmetric energy profile and (2) a po- 
tentially stable (at least in the vicinity of the linear 
spectrum's optical band edge) breather that is cen- 
tered on a light particle and is characterized by an 
asymmetric energy profile. In the present work, we 
focus on the latter, more "robust" solution, and we 
compare it with its counterpart at the interface. It 
should be noted that such comparisons have recently 
been a focal point of interest in other fields, including 
e.g. nonlinear optics where such modes, the so-called 
nonlinear Tamm states, were proposed theoretically 
[28, 29], but also observed experimentally [30-f32jl in 
both one- and higher-dimensional settings. 



FIG. 1: The schematic set-up for the bead 
configuration: steel - aluminum - . . . steel - 
aluminum - brass - brass - ... - brass 



More specifically, we investigate the existence, sta- 
bility, and dynamics of discrete breathers in a com- 
pressed granular chain consisting of a diatomic part 
and a monoatomic part (cf. also with the experimen- 
tal context of Ref. 



where an interface between 
two monoatomic chains was considered). We are in- 
terested in such localized excitations at the interface 
between the two parts of the chain. We first present 
the theoretical model, then analyze its linear spec- 
trum; using the latter as a starting point, and per- 
forming a continuation of the frequency parameter 
within the linear spectral gap, we are able to iden- 
tify a family of nonlinear states localized at the in- 
terface. We present the stability of these states and 
compare/contrast it with that of the bulk structures 
in a diatomic chain identified in [lH , and 22 [ . We 
remarkably find that these interface breathers are 
more robust than their diatomic bulk analogs. Next, 
we present another family of interface breathers that 
are shifted away from the interface by two lattice 
sites (a unit cell in the diatomic chain), and cease to 
exist as their frequency approaches the optical band 
edge. Lastly, we study the dynamics of this model, 
investigating the minimum requirements for gener- 
ating these breathers in experimentally relevant set- 
tings, as well as examining the role of the interface in 
interactions with highly nonlinear, traveling solitary 
waves that are also supported by the model. 



THEORETICAL SETUP 

We will focus our considerations on a one- 
dimensional array of beads, uniaxially compressed 
by a force Fq, which is modeled by zero-tension non- 
linear spring-mass system. The presence of a pre- 
compressive force Fq , ensures the existence of an un- 
derlying linear model and hence a linear spectrum. 
Note that the chain and the precompressive force 
will be assumed to be horizontal, so that the effects 
of gravity will be absent. The equation of motion 
for the displacement m„ of the nth bead (away from 



the chain boundaries) will be given by [l| 



^n — l,n [^n— l.n 



(u„+l - Un)W, (1) 



where An^n+i is the elastic coefficient between beads 
n and rt -I- 1, and depends on the exponent p and the 
geometric and elastic properties of the beads n and 
n + 1. In the Hertzian setting, the exponent takes 
the value p = 3/2 and the coefficient An n+i is given 

by 



An^n+l — 



(2) 



where r„, £"„, and Vn denote, respectively, the ra- 
dius, elastic (Young's) modulus, and Poisson ratio 
of the nth bead. A„ „+i = (i^o/^n,n+i)^^'^ is the 
amount of static overlap of the beads n and n + 1 
(due to the precompression) under equilibrium con- 
ditions. The bracket [s]+ of Eq. ([1]) takes the value 
s if s > and the value if s < 0. This accounts for 
the fact that if two beads are not in contact, they 
do not exert any force on each other (zero-tension). 
Each bead (away from the boundary) has a local 
energy density associated with it Js^] given by 



= ]^rnnul^ + ]^[Vn-l{Un-Un-l) + Vn{u„+i 



where the interaction potential V^„(u„+i — it,, 
given by 



■u„)], 
(3) 



5/2 



An,n+l^J^j^l{Un+l ~ Un) 



(4) 



For finite chains of length iV, we have implemented 
free boundary conditions thus, the equations of mo- 
tion for the displacements of the end-beads read: 



miili = Fq - Ai,2[Ai,2 - {U2 - Ul)] 



3/2 



(5) 



rriNUN = An_i^n[An-i,n - [un - un-i)]+^ - Fq. 

(6) 

More specifically, we consider a chain of alternat- 
ing beads (i.e., a diatomic lattice with the beads 
having the same radius but being made of differ- 
ent materials) and a chain of a third type of bead 
(monoatomic) adjacent to the diatomic fraction of 
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Material Steel Aluminum 


Brass 


2.5^ 


Radius (mm) 4.76 


4.76 


4.76 


2 


Density (kg/m^) 8027.17 


2700 


8500 


^1.5 

N 

X 


Elastic modulus (Pa) 1.93 x lO" 


7 X lO^" 


7.6 X 10"' 


Poisson ratio 0.3 


0.35 


0.414 


" 1 


TABLE I: The material properties used 


in our 


0.5 


simulations. 
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the full lattice. Figure [T] shows the schematic of the 
bead configuration. In the first half of the chain, 
the alternating beads used consist of steel and alu- 
minum; the second half of the chain is composed 
entirely of brass beads. This specific configura- 
tion and these materials were chosen based on their 
straightforward availability for experimental realiza- 
tions, and equally (or more) importantly, because 
of the properties of the underlying linear spectrum 
which enable the existence of localized vibrational 
mode at the interface between the diatomic and the 
monoatomic part of the chain. We will discuss these 
linear properties in more detail in a later section. A 
summary of the material elastic and geometric char- 
acteristics used is given in Table HI 



RESULTS 
Numerical methods 

The numerical integration of the system was per- 
formed using the ode 113 function in Matlab, which 
is a variable order Adams-Bashforth-Moulton solver, 
and allows for stringent error control. To calculate 
the numerically exact (up to a prescribed relative 
tolerance, usually 10~^°) profile of each DB solu- 
tion, we used a Newton-Raphson solver to calculate 
the fixed-points of a Poincare map with a frequency 
/ and, thus, a period T = 1//, 

i?(w(0)) EE u(0) - u{T) = 0. 

Where, we add the additional constraints that m„ = 
in order to "pin" the solution to a point on 
the Poincare map. The Jacobian for the Newton's 
method is given by / — V{T), where V{T) is the 
monodromy matrix, which we obtain by integrating 





Eigenfrequency Index n 

FIG. 2: (Color online) In the left panel, the 
eigenfrequencies of the configuration shown in Fig. 
[T]are shown for 80 beads. The right four panels 
show the first, second, third, and fourth optical 
eigenmodes. For clarity, steel beads are shown as 
gray circles, aluminum beads are shown as blue 
diamonds, and brass beads are shown as brown 
squares. 



the variational equation 

simultaneously with the equations of motion (fTj). 
where J is the Jacobian of these nonlinear equations, 
and V{0) = I. See [s^, and the references therein 
for more details. The original initial guess for such 
an interface DB consisted of the eigenmode of the 
lower cutoff edge of the optical band of the linear 
spectrum. However, in the case of the bulk DBs, 
there was no eigenmode that was localized only in 
the bulk, and we needed to first obtain one solution 
from a purely diatomic chain, which we then embed- 
ded in the full chain as our initial guess. We then 
performed a continuation of the frequency through 
the linear spectrum band gap 37[. This was done 
both for the interface and for the bulk DBs. It 
should additionally be noted that, upon convergence 
of the relevant iteration, the eigenvalues of the mon- 
odromy matrix provide us with the Floquet multi- 
pliers determining the linear stability of the corre- 
sponding solution. The presence of eigenvalues of 
unit modulus (i.e., of multipliers that reside on the 
unit circle) is a signature of linear stability of the 
DB modes, while the existence of eigenvalues with 
modulus greater than one is a signature of (linear) 
instability, as discussed also in more detail below. 
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Underlying linear spectrum 

The eigenfrequencies of the configuration of Fig. 
[U are shown in Fig. [5] These are computed by 
eigenanalysis of the linearized Eqs. ([1]), ([5]), and 
(|6|). The linear spectrum of the diatomic part of 
the chain contains two branches {acoustic and opti- 
cal) and possesses a gap between the upper cutoff 
frequency uif = \j2KsAl'fns of the acoustic branch 
and the lower cutoff frequency uj^ ~ yj2KsA/'mA 
of the optical branch, where ms (jua) is the mass 
of the steel (aluminum) beads. The linear stiffness 

o 2 /3 1 /3 

of the steel-aluminum contact is Ksa — ^^sa^o ■ 
The upper cutoff frequency of the optical band is 
located at a;| = y/2KsA{^/mA + l/mg). The lin- 
ear spectrum of the monoatomic part contains only 
an acoustic branch with frequencies from to = 
y^AKBs/iTT-B: whcrc niB is the mass of the brass 
beads and Kbb — ^^^bb^o^^ linear stiffness of 
the brass-brass contact. 

For the material properties mentioned above and 
for a precompressive force of Fq = 22.4 N, the cut- 
off frequencies obtain the following values: /™ = 
15.21kHz, = 11.96kHz, = 20.6kHz, and 
/3 = 23.8kHz. Since the upper edge of the acoustic 
band of the monoatomic part of the chain is less than 
the lower edge of the optical band for the diatomic 
part of the chain, we obtain a gap in the spectrum 
of the union of the two parts of the chain, between 
/™ = 15.21kHz and = 20.6kHz. The materials 
we used were chosen in such a way that the acoustic 
band of the monoatomic part of the chain did not 
completely overshadow the band gap in the diatomic 
part of the chain. We can widen or narrow the gap 
based on the materials chosen and on the magnitude 
of the precompressive force. 

The eigenmode from the bottom edge of the opti- 
cal band was used as an initial guess in our Newton's 
method solver. An example guess, for 80 beads, is 
shown in the top of the right panels of Fig. [2] Later, 
the second and third optical eigenmodes were used 
as an initial guess in continuations through the gap. 
These eigenmodes are shown in Fig. [2] as well. 

As one can see, the eigenmode from the bottom 
edge of the optical band is localized right on the in- 
terface between the two halves of the chain. This 
localization is caused by the use of free boundary 
condition (BC) on the left and the interface between 
an aluminum and brass bead. In an infinite di- 
atomic chain, the lower optical mode has a profile 
in which the heavy particles (steel in our case) are 
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FIG. 3: (Color online) Examples of four solutions. 
The left two are located at the interface, the right 
two are in the bulk; the top two are closer to the 
optical band edge, while the bottom two are deeper 

within the spectral gap, and thus exhibit a more 
pronounced localization. For clarity, steel beads are 
shown as gray circles, aluminum beads are shown 
as blue diamonds, and brass beads are shown as 
brown squares. 



at rest while the light particles (aluminum) move in 
an alternating pattern with the same amplitude. In 
a semi-infinite diatomic chain, with fixed left BC, 
[u{0) = 0] when the first particle is light, the form of 
this mode is exactly the same as that of the inifinite 
chain since the "wall", particle 0, is a heavy parti- 
cle at rest as it should be. However, using a heavy 
particle at the first site, or alternatively, a free left 
BC [u(0) ~ u{l)], leads to a lower optical mode in 
which the light masses move again in an alternate 
pattern but now with a decreasing amplitude as we 
approach the left part of the chain. In our configu- 
ration, the diatomic part of the chain has a left free 
BC and ends with a light particle at the interface. 
The latter acts more like a wall, rather than a free 
end; the brass bead is relatively heavy, with a mass 
close to that of a steel bead, and during the exci- 
tation of the lower optical mode, the displacement 
of the brass bead is very small around the equilib- 
rium since this frequency belongs to the gap of the 
underlying monoatomic spectrum of the brass. On 
the other hand, according to the above explanation, 
an interface between a steel and brass bead leads to 
a mode with substantially decayed amplitude at the 
interface. Thus, a configuration with that kind of 
interface cannot lead to an interface breather that 
bifurcates from the linear limit of the lower edge of 
the optical branch. 
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FIG. 4: Top: the magnitudes of the Floquet 
muhiphers for the bulk (left) and mterface (right) 
breathers are shown as a function of frequency. In 
both cases, the real muhiphers are shown in black, 
while the complex multipliers are shown in gray. 

Bottom: the energies of the bulk (left) and 
interface (right) breathers are shown as a function 
of frequency. Note that the vertical scales are the 
same in each row, and the horizontal scales are the 
same in each column. 



Existence and stability of interface discrete 
breathers 



Using the eigenmode from the lower edge of the 
optical band as an initial guess, we performed a con- 
tinuation as a function of the frequency through the 
linear spectral gap. This resulted in the identifica- 
tion of the DBs that were located at the interface 
between the diatomic and monoatomic chains, to 
which we will refer as interface DBs. In addition, for 
comparison, we used the corresponding linear eigen- 
mode of a purely diatomic chain to solve for DBs in 
the bulk of the diatomic part of the chain (i.e. away 
from the interface and boundaries) . We will refer to 
the latter class of solutions as bulk DBs. Some exam- 
ples of typical obtained DB profiles are shown in Fig. 
[H The further into the middle of the gap the fre- 
quency of the breather is located, the more localized 
the relevant structure becomes (a feature shared by 
the interface and bulk waveforms). Approaching the 
acoustic branch, both interface and bulk waveforms 
delocalize and obtain non-zero oscillating tails due 
to resonance with the upper acoustic cutoff mode 
of the monoatomic brass part. More details about 



this kind of delocalization for the case of the bulk 
breathers in a diatomic chain can be found in [2l|. 

On the other hand, the linear stability analysis 
illustrated one major difference between the bulk 
breathers and their interface analogs. For both types 
of DBs, the total energy of each breather solution 
(E = ^») computed as a function of the 

corresponding frequency (during the frequency con- 
tinuation). In the plot of energy vs. frequency for 
the bulk case, there is a change in the monotonicity. 
This has been found to correspond to a change in 
the stability of the solutions (shown in Fig. U) (no- 
tice that this type of result suggests more broadly 
that a theorem of this type could be mathemati- 
cally proved, by analogy to the famous Vakhitov- 
Kolokolov criterion for stability of solutions in non- 
linear Schrodinger type models [s^). More specif- 
ically, when the energy of the solution is decreas- 
ing with respect to its frequency, the corresponding 
breather is either linearly stable (e.g., very close to 
the optical band edge) or only oscillatorily unstable. 
The latter oscillatory instabilities are due to quar- 
tets of Floquet multipliers, with non-zero imaginary 
parts, that depart from the unit circle. Based on our 
observations, the magnitude of oscillatory instabil- 
ities decreases as the length of the chain increases, 
so these are essentially weak instabilities that only 
manifest themselves over very long time scales (much 
longer than the potential experimental lifetime of 
the breather [22|). One important note about these 
bulk DBs is that, unlike the case of a purely diatomic 
chain, they do not bifurcate from the optical band, 
and because of this, we could not perform the contin- 
uation all the way to the top of the gap (once again, 
see [2l| for more details about the purely diatomic 
case). On the other hand, when the energy of the 
solution is increasing with respect to the frequency, 
the corresponding breather shows a definitive signa- 
ture of a real instability, corresponding to a pair of 
Floquet multipliers bifurcating off of (1,0) and on 
to the real axis. See Fig. |4]for a more detailed look 
at the Floquet multipliers. Remarkably, in the plot 
of energy vs. frequency for the interface breathers 
(shown in Fig. |4]), there is no such change in mono- 
tonicity, leading to the conclusion that throughout 
the parametric continuation, the interface breathers 
are either stable or merely subject to weak oscilla- 
tory instabilities. This suggests that the interface is 
acting as a stabilizer for the breathers and may be 
useful for experimental purposes of ensuring the ro- 
bustness/observability of such structures especially 
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FIG. 5: (Color online) Top: the magnitude of the 
Floquet multipliers for the double (left) and triple 

(right) breathers is shown as a function of 
frequency. In both cases, the real multipliers are 
shown in black, while the complex multipliers are 
shown in gray. Bottom: example profiles of the 
double (left) and triple (right) breathers for a 
frequency of 18.3 kHz. For clarity, steel beads are 
shown as gray circles, aluminum beads are shown 
as blue diamonds, and brass beads are shown as 
brown squares. 



in the vicinity of interfaces. It is even more remark- 
able that similar (yet reverse in their nature) stabil- 
ity modification conclusions have been reported in 
the literature of nonlinear Schrodinger type systems, 
see e.g., [11], and [s^l, where it was found that sur- 
face breathers may be less robust (and subject to an 
interval of instability) than their bulk counterparts. 
However, it should be noted that the latter find- 
ings are for lattices that even in their monoatomic 
realization possess localized solutions of this type 
(whereas for the granular chains considered here, a 
diatomic setting is necessary for the existence of bulk 
breathers within the linear spectral gap 22, 2^). 



Continuation of second and third optical 
eigenmode 

The first optical eigenmode was not the only 
eigenmode that may lead to localization upon con- 
tinuation within the spectral gap. Figure [2] shows 
that the second and third eigenmodes also led, upon 
corresponding continuations, to states with two or 



three, respectively, breather structures, with one lo- 
calized at the interface in each case. This trend 
can be extended to higher eigenmodes which simi- 
larly lead to corresponding multi-breather structures 
within the spectral gap. We performed a continua- 
tion through the gap using the second, and then the 
third, optical eigenmode as an initial guess (once 
again, with 80 beads). Figure [5] shows a sample pro- 
file from each of these continuations (as well as their 
parametric stability properties). The resulting so- 
lutions contained two, and three, breathers, respec- 
tively, which can be clearly discerned to develop near 
the regions of maximal displacement in the respec- 
tive eigenmodes. Once again, the linear stability of 
each solution was analyzed using the Floquet multi- 
pliers. A detailed view of the obtained multipliers is 
shown for the double and triple breathers in Fig. [S] 
These double and triple breathers have intervals of 
instability similar to the bulk breathers, which im- 
plies that it is the bulk breathers "within the struc- 
ture" that become unstable, while the localization at 
the interface may be expected to persist. This is con- 
firmed by dynamical simulations of these breathers, 
in which we make the following observations. When 
we perturb an unstable double, or triple breather, 
the unstable bulk part radiates some energy and de- 
cays to one of the stable bulk breathers (with a lower 
frequency) , while the interface breather remains un- 
changed. Then the two (or three) breathers oscillate 
out of sync, until oscillatory instabilities disperse the 
breathers: usually the bulk breathers first, followed 
by the interface breather. 



Another family of discrete breathers 

Although in the above considerations, we exam- 
ined the different families of breathers arising from 
the linear limit, an interesting question concerns 
whether additional families of surface breathers re- 
sembling the ones above may exist, shifted away 
from the interface; such families may, in fact, cease 
to exist before the linear limit is reached since they 
do not degenerate into linear modes. For example, 
as we will see in the dynamical evolution simula- 
tions below, under certain conditions (of interac- 
tion of a traveling compression pulse with an in- 
terface breather), a different type of surface mode 
may be excited. In this mode, the maximal dis- 
placement site is shifted away from the interface by 
one unit cell (two lattice sites). Using an appropri- 
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FIG. 6: (Color online) Results for a continuation of 
a DB located one unit cell from the interface. Top: 
the magnitude of the Floquet multipliers for each 
DB is shown as a function of frequency. The purely 
real multipliers are highlighted in black. Middle: 
the energy of each DB is shown as a function of 
frequency. Bottom: a sample DB solution from this 
family is shown. Note that for visual clarity, brass 
is shown in brown (squares), steel is shown in gray 
(circles), and aluminum is shown in blue 
(diamonds). 



ate initial guess of a breather with an appropriately 
shifted center, we performed a continuation through- 
out the frequency gap, and obtained such a family 
of breathers which were localized one unit cell away 
from the interface. The bottom panel in Fig. IH] 
shows a sample profile from this continuation, while 
the top two panels show the magnitude of the Flo- 
quet multipliers, and the energy as a function of the 
DB frequency. It is interesting to note that this fam- 
ily of breathers does not bifurcate from the optical 
band, as the other family of interface DBs does. In 
fact, as the frequency of these DBs approaches the 
optical band, they develop a real instability (cf. also 
the associated change of monotonicity in the energy- 
frequency plot in Fig. [S]). These real instabilities 
increase rapidly until the family undergoes a saddle- 




FIG. 7: (Color online) Top left: a spacetime plot of 
the local energy density for a compression pulse 
passing through the interface. There is a barely 
visible reflection amounting to about 2% of the 
total energy. The lower left panel shows a different 

view of the same simulation. Here, each line 
represents the total energy in a certain part of the 
chain. The dashed line (blue) represents the first 
half of the chain, the solid line (red) represents the 
interface region, and the dash-dotted line (green) 
represents the second half of the chain. There is no 
significant amount of energy trapped in the 
interface. The right two panels show a similar 

picture, except that in this case, there is an 
interface breather for the pulse to interact with. 
The pulse causes the breather to shift by 1 cell. 
Note that the time-scale is the same for the top 
and bottom panels, and the color scale for the top 
panels is truncated at 1/20 the maximum, in order 
for the refiections to be visible. 



center bifurcation, and ceases to exist. 



Discrete breather and nonlinear interface 
dynamics 

We now turn to direct numerical simulations in or- 
der to examine the nonlinear dynamics of the chain 
(and, in particular, of the interface), as well as to 
assess the role of the discrete interface breathers 
therein. 

Our first set of numerical experiments involved 
scattering events. The aim of these was to inves- 
tigate the behavior of a compression pulse as it in- 
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FIG. 8: (Color online) Top: the destruction of a 
bulk breather due to interaction with a compression 
pulse. Bottom: the destruction of an interface 
breather due to the oscillatory instabilities. 
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FIG. 9: (Color online) The frequency of the 
induced breather (or linear mode) when the middle 
site (to the left of the interface) was excited with a 
negative displacement is shown as a function of the 
magnitude of the initial excitation. Recall that A 
is the amount of overlap between the beads at 
equilibrium. The dashed lines represent the 
acoustic and optical band edges. 



teracted with the interface, in the absence, as well as 
in the presence of an excited breather. We excited 
the first bead with a striker to send a compression 
pulse through the chain. Many of these simulations 
were conducted over a wide range of velocities for 
the striker, and Fig. [7] shows typical examples within 
this class of numerical simulations. In this particu- 
lar simulation, the velocity of the striker was about 
0.94 m/s, and produced a solitary wave with an am- 
plitude of about 2.5 x 10~^ m and a maximum dy- 
namical force of approximately 138 N, but similar re- 
sults were obtained for lower initial striker velocities. 
Note that in this figure, the space-time contour plots 
of the local energy density (defined above), rather 
than the displacements, are shown. The leftmost 
plots in Fig. [7] show a compression pulse colliding 
with the interface; although we observe some mild 
reflection (about 2%; see also below) at the interface, 
it is generally found that the pulse does not excite a 
localized mode at the interface. The rightmost plots 
in Fig. [7] show a compression pulse colliding with a 
pre-existing breather at the interface. Surprisingly, 
there does not appear to be a significant difference in 
the behavior of the pulse in this second case, but it 
does illustrate that under certain circumstances, the 
breather actually shifts away from the interface upon 
collision with the pulse. The dynamics here clearly 
appears to lead to a member of the family of surface 
discrete breathers studied in the previous section. 
After further tests, we conclude that whether or not 
a breather shifts away from the interface appears to 



depend on its linear stability. It seems to be the case 
that a compression pulse of the right magnitude will 
shift an interface breather away from the interface 
only if the breather has a pair of Floquet multipliers 
that split off the unit circle at the (—1, 0) point of the 
unit circle. I.e., such subharmonic instabilities can 
lead to potential mobility of the interface breathers 
under the action of suitably strong perturbations. In 
order to measure the transmission of energy past the 
interface, we measured the change in energy in each 
part of the chain (the first diatomic half, the inter- 
face, and the monoatomic second half), and divided 
by the total energy. The bottom two panels of Fig. 
[7] show a typical example of how the energy propa- 
gated through the chain (in this case, for a chain of 
1000 beads). The left hand part is the total energy 
in beads 1-497, the interface part is the total energy 
in beads 498-502, and the right hand part is the total 
energy in beads 503-1000. The reflection of the soli- 
tary wave from the interface is observed as the small 
dip in the dashed line after the pulse passes through 
the interface and only amounts to about 2% of the 
total energy. Notice that although the reflection is 
fairly weak, the fraction of energy remaining on the 
left part of the domain is more substantial; this is 
due to the compression pulse formation process from 
the striker, which has led to residual "radiation" of 
phonons in its wake (not shown here, since for clarity 
we only show the part of the domain in the vicinity 
of the interface). 

On the other hand, we note in passing that, in 
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the case of the unstable bulk breathers (through a 
real pair of Floquet multipliers), their interaction 
with a compression pulse has been observed to lead 
to the ultimate destruction (i.e., dispersion) of the 
breather mode, which is shown in Fig. [Sj Lastly, 
even in the absence of the solitary wave, the inter- 
face breathers will eventually decay, due to the os- 
cillatory instabilities, as shown in Fig. [S] This is 
a similar phenomenon to what is observed also for 
bulk breathers when they are subject to the same 
type of instabilities [2lj . 

Another set of performed numerical experiments 
concerned the realization of the minimum initial 
excitation needed to produce an interface breather 
with a frequency in the spectral gap. From the sim- 
ulations above, it was clear that compression pulses 
did not excite such a breather, regardless of which 
end they started from. This is similar to what is 
shown in the bottom left panel of Fig. [7] where it 
shows the interface region has a brief spike as the 
pulse passes through the interface, but then returns 
to zero. Exciting the middle site with a positive dis- 
placement was also unsuccessful. The positive dis- 
placement excites the monoatomic part of the chain, 
namely acoustic modes up to 15.21kHz. But this 
part of the chain is dispersive and the nonlinear- 
ity can not support localized vibrations. Thus, in 
all cases of this type, the energy was dispersed into 
linear waves and nonlinear solitary waves, and no 
substantial fraction thereof was found to be trapped 
at the interface. 

On the other hand, when the middle site was ex- 
cited with a negative displacement of sufficient am- 
plitude, we observed the formation of discrete inter- 
face breathers. The larger the amplitude of the ini- 
tial excitation, the lower the frequency into the gap 
the resulting breather is found to be (Fig. IH] shows 
the frequency of the resulting breather as a function 
of the initial excitation). The top panels of Fig. [TOl 
show the results of a simulation in which the central 
site was excited to about 1.5 A. It can be observed 
that, after some initial solitary waves are released, 
a robust interface breather is formed almost imme- 
diately. The resulting breather, which, in this case, 
has a profile shown in Fig. IIOI has a frequency of 
about 16.193kHz, and matches almost perfectly the 
exact breather of that frequency identified through 
the Newton's method. These observations may be of 
use in trying to produce these modes experimentally. 
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FIG. 10: (Color online) Top left: The space-time 
graph showing the energy density for a chain of 
1000 beads in which the 500th bead was excited to 
about 1.5 A. After some initial radiation, a DB is 
excited at the interface. Top right: A zoom of the 
top left panel, showing that there is indeed a DB 
located at the interface. Bottom: the profile of the 
excited breather is shown in conjunction with the 
exact breather for the corresponding frequency. 



CONCLUSIONS AND FUTURE 
CHALLENGES 



In the present work, we considered the possibil- 
ity of formation of discrete breathers at the inter- 
face between a diatomic chain and a monoatomic 
chain (of same radii beads but made of different 
materials). We chose this particular configuration 
due to its underlying linear properties and eluci- 
dated how (and under what conditions) these give 
rise to a linear spectral gap whose optical band cut- 
off mode enables a continuation towards a discrete 
breather state. The continuation of such modes was 
performed throughout the spectral gap and the ap- 
parent robustness of these structures was elucidated. 
It was, in fact, illustrated that these modes are ar- 
guably more robust than their bulk counterparts 
possessing no parametric intervals of non-unit Flo- 
quet multipliers along the real axis (as their bulk 
counterparts do), but being only subject to weak 
oscillatory instabilities for finite-size chains. This 
should imply that similarly to what was shown for 
their bulk counterparts in [2^ . such modes should 
be possible to experimentally produce in granular 
chains. In this general direction of interfaces, the 
work of [13] has experimentally considered the prop- 
agation of waves at the interface of two different 
granular chains. More recently, in the same spirit. 
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the work of 27| has focused on the dynamics of com- 
pression pulses at the interface of an uncompressed 
monoatomic chain with a rod. In our setting, in 
addition to illustrating the existence and stability 
of breather interface modes, we also examined how 
they behave dynamically and how they interact with 
compression pulses. While the latter are unable to 
excite such interface breathers, we have illustrated 
under what conditions (and negative displacements 
of the interface bead) such modes may arise. Fur- 
thermore, we have shown that their interaction with 
compression pulses may lead to their shift from the 
interface. On the other hand, it is interesting to 
mention that we have also revealed multi-breather 
states containing both breathers located at the in- 
terface, as well as (one or more) bulk ones. These 
stem from higher order linear modes of the spec- 
tral problem involving the interface. Lastly, we have 
shown the existence of a family of interface DBs, 
which does not bifurcate from the optical band, as 
do the other families of interface DBs. 

Our results suggest a number of interesting direc- 
tions for future explorations. On the one hand, they 
raise the question of appropriate quantitative diag- 
nostics for the examination of interaction of pulses 
(or breathers) with interfaces. On the other hand, a 
natural topic that emerges concerns the conditions 
under which a propagating excitation (e.g. a nearly 
compact pulse) may produce energy trapping at an 
interface similar to the ones used herein. Finally, 
it would be interesting to examine similar phenom- 
ena (and interfacial dynamics) in other classes of 
models, such as a binary interfacing with a unary 
waveguide array modeled by the discrete nonlinear 
Schrodinger equation to examine the potential gen- 
erality of the conclusions obtained herein [in particu- 
lar, ones about the enhanced stability of the surface 
modes]. Such investigations are presently underway 
and will be reported in future publications. 
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